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NOTICE 


m i n ar y an  a I ys  i •=-  of 
i ent  man i tar i ng  data 

The  following  report  is  based  on  a preliminary  analysis  of 
data  collected  during  the  first  two  years  of  a five-year  data 
collection  program.  The  conclusions  of  this  report  are  therefore 
subject  to  change  after  completion  of  the  entire  five-year 
program  and  a thorough  analysis  of  all  data. 
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INTRODUCTION 


The  1987  Chesapeake  Bay  Agreement  states  that  "the 
improvement  and  maintenance  of  water  quality  is  the  single  most 
critical  element  in  the  overall  restoration  and  protection  of  the 
bay."  Improvement  in  water  quality  includes  "the  reduction  of 
nutrients  from  both  point  and  nonpoint  sources."  Agreement  was 
reached  to  set  a goal  for  achieving  a 40%  reduction  in  nutrients 
entering  the  bay  system  by  the  year  2000.  Data  that  is  useful  in 
developing  a plan  to  achieve  that  goal  has  been  and  is  being 
collected  for  the  Chesapeake  Bay  Program  by  a joint  program 
directed  by  the  Susquehanna  River  Basin  Commission. 

The  objective  of  this  joint  Susquehanna  River  Basin 
Commission  (SRBC) /Pennsylvania  Department  of  Environmental 
Resources  (PaDER),  Bureau  of  Soil  & Water  Conservation  and  Bureau 
of  Laboratories/U . S . Geological  Survey  (USGS)  program  is  to 
collect,  collate,  analyze  and  disseminate  nutrient  data  for  the 
Susquehanna  River  and  selected  tributaries.  This  program 
establishes  a sound  database  which  can  be  used  to  effectively 
plan  and  implement  both  immediate  and  long-range  nutrient 
reduction  efforts  and  support  modeling  activities. 

These  data  were  also  used  to  estimate  the  nutrient  loads  for 
this  report. 
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SAMPLE  COLLECTION  AND  ANALYSIS 


Sample  collection  began  at  twelve  sites  in  the  Susquehanna 
River  basin  in  October,  1984.  A listing  of  the  sampling  sites 
and  their  location  is  on  the  following  page.  Yearly  sampling 
consists  of  collecting  twelve  monthly  baseflow  samples  as  well  as 
storm  samples  during  each  of  the  four  seasons.  Storm  coverage 
includes  sample  collection  on  both  the  rising  and  falling  stream 
stage . 

The  depth-integrated,  multivertical,  composited  samples  are 
brought  to  the  PaDER  lab  usually  within  24  hours  of  collection. 
Nitrite  plus  nitrate  -N  and  ammonia  -N  are  normally  analyzed 
within  24  hours  of  receipt  of  samples.  Only  EPA  approved  methods 
of  chemical  analyses  are  used  by  the  PaDER  lab  for  the  programs 
samples . 

Chemical  constituents  analyzed  include: 
nitrite  + nitrate  - N,  total 
ammonia  - N,  total 
Kjeldahl  N,  total  and  dissolved 
Phosphorus,  total  and  dissolved 
Orthophosphate  - P,  dissolved 
Organic  carbon,  total 

Concentrations  of  suspended  sediment  are  determined  for 
all  samples  at  the  USGS  sediment  laboratory. 


2 


SECTION  OF 

SUSQUEHANNA  RIVER  BASIN  SHOWING  LOCATIONS 
OF  NUTRIENT  MONITORING  SITES 
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LEGEND 

© EXISTING  U.S.G.S.  GAGE 
A NEW  GAGE  LOCATION  / SAMPLING  SITE 


The  data  collected  under  this  program  were  used  to  estimate 
the  nutrient  loads  for  the  entire  Susquehanna  River  Basin  as  well 
as  at  the  12  sites  within  the  basin  for  Calendar  Years  1985  and 
1986  . 


The  annual  nutrient  load  for  the  Susquehanna  River  Basin  are 
shown  in  the  table  below. 


ANNUAL  NUTRIENT  LOADS  FOR  THE  SUSQUEHANNA  RIVER  BASIN 


CY-1985 
(Tons ) 

CY-1986 
(Tons ) 

Organic  Nitrogen 

17,215 

28,548 

Inorganic  Nitrogen 

47,190 

55,740 

Total  Nitrogen 

64,405 

84,827 

Total  Phosphorus 

1,992* 

4,598 

Suspended  Sediment 

776,017* 

3,626,683 

Dissolved  Phosphorus_ 

1,197 

1,711 

Total  Nitrite  plus  Nitrate  Nitrogen 

45,123 

52,825 

Total  Ammonia  as  Nitrogen 

2,066 

2,915 

Dissolved  Kjeldahl  Nitrogen 

13,742 

20,298 

Total  Kjeldahl  Nitrogen 

20,611 

33,210 

Total  Organic  Carbon 

103,043 

162,800 

* Results  questionable. 


Annual  nutrient  loads  fluctuated  from  CY-1985  to  CY-1986. 
CY-1985  annual  loads  were  computed  using  primarily  baseflow  data 
while  CY-1986  data  included  a number  of  storms.  This  mirrors  the 
difference  in  runoff  between  the  two  years.  The  lower  runoff  for 
1985  is  reflected  by  the  fact  that  the  average  discharge  per 
square  mile  was  generally  about  1 cfsm  while,  except  for  Codorus 
Creek,  the  average  in  1986  was  nearly  1.5  cfsm.  Peak  discharges 
were  also  generally  lower  in  1985  than  in  1986  . Discharge  for 
CY-1986  was  slightly  above  average  based  on  data  for  the  USGS 
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gage  at  Marrietta,  PA.  The  average  discharge  for  1986  was  38,280 
cfs.  compared  to  the  55-year  average  of  36,900  cfs.  Also, 
Williams  and  Reed  (1972)  reported  that  under  average  runoff 
conditions,  the  Susquehanna  River  transports  about  3 million  tons 
of  sediment  annually.  Approximately  3.6  million  tons  of  sediment 
were  estimated  to  have  been  transported  in  1986. 


PERCENTAGE  OF  BASIN  ANNUAL  NUTRIENT  LOAD 
ORIGINATING  UPSTREAM  OF  HARRISBURG 


CY-1985 

CY-1986 

Organic  Nitrogen 

79 

83 

Inorganic  Nitrogen 

83 

77 

Total  Nitrogen 

82 

78 

Total  Phosphorus 

42 

78 

Suspended  Sediment 

41 

86 

Dissolved  Phosphorus 

59 

69 

Total  Nitrite  plus  Nitrate  Nitrogen 

83 

77 

Total  Ammonia  as  Nitrogen 

68 

70 

Dissolved  Kjeldahl  Nitrogen 

82 

82 

Total  Kjeldahl  Nitrogen 

73 

78 

Total  Organic  Carbon 

78 

80 

Drainage  area  above  Harrisburg  = 88%  of  entire  Susquehanna 
River  basin 


With  the  exception  of  suspended  sediment  and  phosphorus  for 
the  CY-1985 , above  table  shows  the  the  percentage  of  annual  loads 
originating  upstream  of  Harrisburg  corresponded  fairly  well  with 
the  percentage  of  drainage  area  upstream  of  Harrisburg.  It  was 
previously  pointed  out  that  flows  in  CY-1985  were  low.  Low  flow, 
particularly  on  the  Main  Stem,  results  in  low  stream  velocities, 
velocities  probably  not  high  enough  to  transport  sediment.  This 
probably  resulted  in  channel  storage  of  sediment  in  the  Main  Stem 
above  Harrisburg.  This  would  also  be  reflected  in  phosphorus 
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which  has  a high  affinity  for  sediment.  This  stored  material 
will  be  available  to  be  scoured  during  higher  flow  years  causing 
much  higher  loadings  than  that  being  transported  by  overland 
runoff. 


PERCENTAGE  OF  ANNUAL  NUTRIENT  LOAD  THAT 
ORIGINATES  UPSTREAM  OF  DANVILLE  & LEWISBURG 


CY-1985 

Danville  | Lewisburg 

Organic  Nitrogen 

68 

32 

Inorganic  Nitrogen 

63 

37 

Total  Nitrogen 

65 

35 

Total  Phosphorus 

75 

25 

Suspended  Sediment 

71 

29 

Dissolved  Phosphorus 

73 

27 

Total  Nitrite  plus  Nitrate  Nitrogen 

66 

34 

Total  Ammonia  as  Nitrogen 

67 

33 

Dissolved  Kjeldahl  Nitrogen 

70 

30 

Total  Kjeldahl  Nitrogen 

64 

36 

Total  Organic  Carbon 

68 

32 

CY-1986 

Danville 

Lewisburg 

Organic  Nitrogen 

67 

33 

Inorganic  Nitrogen 

63 

37 

Total  Nitrogen 

65 

35 

Total  Phosphorus 

75 

25 

Suspended  Sediment 

79 

21 

Dissolved  Phosphorus 

63 

37 

Total  Nitrite  plus  Nitrate  Nitrogen 

63 

37 

Total  Ammonia  as  Nitrogen 

70 

30 

Dissolved  Kjeldahl  Nitrogen 

68 

32 

Total  Kjeldahl  Nitrogen 

66 

34 

Total  Organic  Carbon 

74 

26 

Percentage  Drainage  Areas  = Danville  62  Lewisburg  38 
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The  data  in  the  above  table  indicate  that  the  loads,  based 
on  being  proportional  to  the  drainage  area,  correspond  very  well. 
In  most  instances,  except  for  the  underlined  data  for  phosphorus, 
suspended  sediment,  and  TOC,  the  loads  were  generally  within  -10% 
of  the  drainage  area.  It  appears  that  in  both  flow  regimes 
phosphorus  was  disproportionately  higher  in  the  Main  Stem 
Susquehanna  River  than  in  the  W.  Branch  Susquehanna  River. 

An  effort  has  been  made  to  establish  that  CY-1985  was  a 
relatively  dry  year  and,  therefore,  should  be  a low  load 
producing  year.  An  anomaly  to  this  perception  is  Codorus  Creek. 
The  following  table  is  a compilation  of  load  data  for  Codorus 
Creek  at  Pleasureville . 


ANNUAL  LOAD  IN  TONS 


YEAR 

ON 

IN 

TN 

TP 

SS 

CFSM 

PEAK  CFS 

1985 

340 

680 

1,020 

100 

55,300 

0.92 

6,600 

1986 

430 

770 

1,200 

70 

33,400 

0.94 

2,030 

As  seen  from  the  table,  inorganic  and  total  nitrogen  (IN,  TN) 
values  are  higher  for  1986  than  for  1985,  yet  total  phosphorus 
and  suspended  sediment  (TP,  SS)  are  lower.  This  apparent  anomaly 
can  be  explained  by  the  stream  discharge.  While  the  average 
discharge  was  slightly  lower  in  1985  than  in  1986  , the  peak 
discharge  was  three  times  greater.  The  significance  of  this  high 
peak  flow  in  1985  is  that  in  four  days  it  carried  44,400  tons  or 
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80%  of  the  annual  sediment  load  while  carrying  only  53  tons  or 
10%  of  the  annual  IN  load  during  the  same  time  period. 
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BASIS  FOR  LOAD  COMPUTATIONS 


Constituent  loads  were  estimated  using  techniques  described 
by  Miller  (1951).  The  load  estimates  were  derived  by  combining 
daily  load  data  obtained  from  the  transport  curve  with  streamflow 
data  taken  from  a flow  duration  curve. 

The  transport  curve  was  obtained  by  converting  the 
instantaneous  constituent  concentration  in  mg/L  into  tons  per 
day  or  pounds  per  day  using  the  equation: 

L = C Q F 

where 

L = constituent  load  in  tons  per  day  or  pounds  per  day 
C = constituent  concentration  in  milligrams  per  liter 
Q = instantaneous  stream  discharge  in  cubic  feet  per  second 
F = 0.0027  for  tons  per  day;  5.39  for  pounds  per  day 

Constituent  loads  were  plotted  against  instantaneous  water 
discharge  on  logarithmic  paper  (for  example,  see  next  page). 
Linear  regression  analysis  were  performed  on  the  logarithmic 
transformations  to  determine  the  relation  between  load  and 
streamflow.  Transport  curves  were  fitted  by  the  method  of  least 
squares.  The  T test  was  performed  to  test  the  slope  for 
significant  difference  from  zero  at  the  95  percent  confidence 
interval.  Regression  equations  and  correlation  coefficients  (r) 
are  listed  with  each  transport  curve. 
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CODORUS  CREEK  ff  YORK:  TOTAL  PHOSPHORUS  CY88 
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X-axia:  Discharge  in  cfa 
Y-axls:  Load  in  lbs/day 
-CYB8:  LOO  (Y)-  -0.849  + i.411KLO0(X) 


The  flow  duration  curve  is  a cumulative  frequency  curve  that 
shows  the  percent  of  time  specified  discharges  were  equalled  or 
exceeded  during  a given  period  (Searcy,  1959). 

Fishel  (1984)  calculated  nutrient  loads  carried  by  the 
Susquehanna  River  at  Harrisburg  using  three  different  methods, 
one  of  which  was  used  in  this  report.  Differences  in  results 
were  usually  less  than  5%  between  the  three  methods. 

The  annual  nutrient  loads  for  the  basin  were  obtained  by 
adding  the  loads  together  for  Susquehanna  River  at  Harrisburg, 
Paxton  Creek,  Swatara  Creek,  W.  Conewago  Creek,  Codorus  Creek, 
and  the  Conestoga  River,  plus  the  estimated  loads  for  the  860 
square  miles  below  Harrisburg  that  were  not  monitored. 
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APPENDIX 

COMPUTED  LOADS  FOR 
CALENDAR  YEARS  1985  AND  1986 
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